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SUMMARY 


The relations between several semiempirical fracture analyses (SEFA) and the 
R-curve concept of fracture mechanics are examined and the conditions for equivalence 
between a SEFA and an R-curve are derived. A hypothetical material is employed to 
study the relation analytically. Equivalent R-curves are developed for several real ma- 
terials using data from the literature. For each SEFA there is an equivalent R-curve 
whose magnitude and shape are determined by the SEFA formulation and its empirical 
parameters. If the R-curve is indeed unique, then the various empirical parameters 
cannot be constant, and vice versa. However, for one SEFA the differences are small 
enough that they may be within the range of normal data scatter for real materials. 


INTRODUCTION 

The relations between several semiempirical fracture analyses and the R-curve 
concept of fracture mechanics are examined. These relations may explain why a semi- 
empirical fracture analysis will yield good results with one set of data and poor results 
with another. They may also indicate which analyses deserve further consideration. 

Over the past decade a number of semiempirical fracture analyses have been pre- 
sented (refs. 1 to 5). These analyses all attempt to correlate failure stresses for pre- 
cracked tension specimens with initial crack length over a range of crack lengths. The 
correlations involve the determination of one (refs . 1 to 3) or two (refs . 4 and 5) em- 
pirical parameters from test data. The parameters are treated as material properties 
which are independent of the specimen and crack configurations but which are functions 
of specimen thickness and such variables as heat treatment and test temperature. The 
analyses often provide very good correlations using some data sets but poor correla- 
tions using other data sets. To date these analyses have only been formulated for and 
applied to test specimen configurations. Thus their applicability to the design of com- 
plex structural configurations is uncertain. 

The progressive development of the R-curve concept has been reviewed in refer- 
ence 6. The concept postulates that, for a given material and thickness, there is a 
unique relation between the amount of stable crack growth under rising load and the 
crack-tip stress intensity factor. This relation is called the crack extension resistance 
curve, or R-curve, and represents the response of the material in the vicinity of the 
crack tip to externally imposed loading. If the R-curve is known, both failure load and 
critical crack length can be predicted (as functions of initial crack length) for any spe- 



cimen or structural configuration for which an appropriate stress intensity analysis is 
available. Thus the R-curve concept appears to be a much more useful method than any 
of the semiempirical analyses . 

If the R-curve for a given material and thickness is available, one can calculate 
residual strength (fracture stress as a function of initial crack length) for any test spe- 
cimen configuration. The converse should also be true. That is, if a relation between 
fracture stress and initial crack length is available, one should be able to calculate the 
corresponding R-curve. This observation was the impetus for the present study, which 
was undertaken to test the following hypotheses: 

(1) For each semiempirical fracture analysis (SEFA) there is an equivalent 
R-curve (ERC) whose magnitude and shape are determined by the SEFA formulation and 
its empirical parameters. The ERC is equivalent in that it predicts exactly the same 
relation between fracture stress and initial crack length as the SEFA. 

(2) A SEFA will correlate residual strength data closely if its ERC closely matches 
the actual R-curve of the material in question and will correlate poorly if the match is 
poor. 

(3) It should be possible to predict the critical crack length in terms of the initial 
crack length and the empirical parameters. 

This report first reviews some characteristics of the R-curve concept when ap- 
plied to finite-width specimens. Next the conditions for equivalence between a semi- 
empirical analysis and an R-curve are derived. A hypothetical material is employed 
to study the relation between R-curves and semiempirical analyses. Finally, equivalent 
R-curves are developed for several real materials using data from the literature. 


SYMBOLS 

a length of single-tip crack or half-length of double-tip crack, a Q + A 

B specimen thickness 

C m empirical parameter (ref. 1) 

2 

E’ effective modulus, E for plane stress or E/(l - v ) for plane strain, where E 

is Young’s modulus and v is Poisson’s ratio 

G^ strain energy release rate 

crack extension resistance 

G c fracture toughness, G^ or G^ at instability condition 

K c empirical parameter (ref. 3) 

Kf empirical parameter (ref. 4) 


2 


Krp c empirical parameter (ref. 5) 

empirical parameter (ref. 2) 

Kj opening-mode stress intensity factor 

m empirical parameter (ref. 4) 

n number of crack tips (one or two) 

W specimen width 

Y stress intensity calibration factor, Kj/ct r* 

a sensitivity factor (eq. (1)) 

A effective crack extension (sum of physical crack extension plus a plastic zone 

correction) 

X relative crack length, na/W 

a stress normal to crack 

o u ultimate tensile strength 

cr yield strength 

ys 

co empirical parameter (ref. 5) 

Subscripts: 

c at critical or instability condition 

o initial value (prior to loading) 


R-CURVE CONCEPT 

The R-curve concept is illustrated schematically in figure 1 for an infinite body 
containing a crack whose original length is 2a Q . The strain energy release rate is 
given by 

2 

P O’ TTcL 

c . 

a E' 

and represents the driving force (per unit thickness) tending to cause crack propagation. 
The material's resistance to crack propagation, G^, is a function of crack extension, A. 
In R-curve analysis the subscripts A, R, and c are customarily used to denote applied 
and resisting forces and critical values, respectively. The R-curve is located with its 
origin at a = a Q . As stress normal to the crack is applied and increased to 90 percent 
of the subsequent critical stress in figure 1, the crack must extend only a small distance 


3 



to develop a large resistance. At this point the crack extension resistance equals the 
driving force and the crack is stable. As the stress is increased, progressively larger 
amounts of crack extension are required to resist the crack driving force. Finally, at 
the critical stress a the driving force curve and the R-curve are tangent. Beyond the 
point of tangeney the driving force increases faster with crack length than does the ma- 
terial's resistance. This instability condition represents the failure of the body. The 
point of tangeney defines the fracture toughness G c and the critical crack length 2a c . 
Since the driving force curve for an infinite body is a straight line, it should be apparent 
that both the fracture toughness G c and the amount of crack extension at instability A c 
increase with increasing original crack length 2a Q . If the R-curve exhibits a plateau, 

G and A may asymptotically approach limit values. 

In simple finite bodies and test specimens, the presence of stress-free boundaries 
results in an additional increase in the crack driving force as the crack extends toward 
a boundary. Thus the slope of the driving force curve increases continuously with in- 
creasing crack length. The instability condition for a typical finite-width specimen is 
shown in figure 2 and is determined as follows. For a given specimen configuration and 
loading, the dimensionless stress intensity calibration factor is defined as 


K t 


Y = 


a-yja. 


fen (A) 


where Kj is the opening-mode elastic stress intensity factor, A = na/W is the relative 
crack length, and n is the number of crack tips (one or two). If we define a dimension- 
less sensitivity factor as 


Y dA 


( 1 ) 


then (for constant stress) the crack driving force curve and its slope are 

E»G a = Y 2 cr 2 a 


A 2 9 

E* — = Y o(l + 2or) 

da 


For convenience, the crack extension resistance curve and its slope are written here as 

g(A) = E'G r 
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g'(A) = E' 


dG^ 

dA 


At the instability point, = G^ and dG^/da = dG^/dA (see fig. 2). If g(A) and 
g'(A) are mathematically describable, the instability point is determined by the simul- 
taneous solution of the following equations: 


E ’ G c = Y c<j<V A c> = 


( 2 ) 


= YV(i+ 2a c ) = g'(A c ) (3) 


The coefficients Y 
functions of the relative 
tion is seldom possible, 
can be done as follows. 


and a are usually expressed as trigonometric or polynomial 
crack length A. As a result, a closed-form simultaneous solu- 
and numerical methods must be used to solve for A c . This 
Dividing equation (2) by equation (3) and rearranging terms give 


E' — 
da 


g(A ) a + A 
_ & v c' _ o c 

g'(A c ) 1+2 q! c 


( 4 ) 


If the functions g(A) and g'(A) and the appropriate equation for a are substituted into 
equation (4) then, for prescribed values of a Q and W, A c is the least positive root of 
equation (4). This root can be found by any of several numerical methods. Next, the 
function g(A) is evaluated at A c to calculate g(A c ). Finally, the fracture stress cr c 
is obtained from equation (2) using the appropriate value of Y . In this report the ex- 
pressions 


Y = 



for the center-crack specimen (ref. 7) and 


Y = A.~ l/2 (1 - Af 3/2 (2 + A) (0.8 66 + 4.64 A - 13.32 A 2 + 14.72 A 3 - 5.6 A 4 ) 


for the compact (tension) specimen (ref. 8; A s 0.2) are used except where noted. 

For a crack in an infinite body, both G c and A c increase with increasing initial 
crack length. In a finite body, as the initial crack length is increased from zero, G c 
and A c increase at first. However, due to the fact that dY/dA continually increases 
with A, both G and A reach maximum values which depend on the specimen width 

t/ L- 
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W and the forms of both the driving force curve and the R-curve „ As the initial crack 
length is increased still further, both G and A begin to decrease. This behavior is 
shown schematically in figure 3 where instability curves are shown for a wide range of 
initial crack lengths . The locus of all instability points is shown by the dashed line . 
From figure 3 it is also apparent that there are pairs of initial crack lengths, say 
(na o )i = 0.2W and (na Q ) ^ ~ 0.7W, which will have the same critical crack extension, 
(A c ) 1 2> an< 3 fracture toughness, (G c )^ 2 * From equation (2), the fracture stresses for 
these initial crack lengths are related by 

(°c^2 _ <^1 A a o^l + ^c^l, 2 
V (a o>2 + ^c)l,2 

Thus there is a relation between fracture stresses for short cracks and long cracks 
which is implicit in the R-curve concept, and this relation is a function of the specimen 
type and the shape of the R-curve. 

If an equation for the R-curve is available, residual strength can be calculated for 
any specimen size and configuration using the procedure which follows equation (4). If 
the R-curve equation has a simple form, it is sometimes possible to develop a dimen- 
sionless residual strength curve. For example, suppose the R-curve is given by 

g(A) = AA b 

where b < 1.0. Then for center-crack specimens the residual strengths are as shown 
in figure 4 for two values of the exponent, b = 0. 5 (a parabola) and b = 0. 2. For this 
case the shape of the residual strength curve depends only on the exponent b. Here it 
is obvious that the R-curve concept implies a particular residual strength curve which 
is a function of the specimen type and the R-curve shape. This implication will be use- 
ful in later comparisons of R-curves and semiempirical analyses. 

In concluding this section it should be emphasized that, in this report, A is the ef- 
fective crack extension. It is the sum of the physical crack extension plus an adjust- 
ment to account for the effect of crack tip plasticity. The nature of the plasticity ad- 
justment, although technically significant, has no influence on the analyses and conclu- 
sions in this report. 


EQUIVALENCY ANALYSIS 

In the preceding section it was shown that, if a mathematical formulation of the 
R-curve is available, fracture stress can be determined as a function of original crack 
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length. In this section it will be shown that, if an equation for fracture stress as a 
function of original crack length is available, the equivalent R-curve can be determined. 

To do this we must differentiate equation (2) with respect to A^, and that operation 
requires some prior consideration. For a given specimen type, there are many combi- 
nations of specimen width and initial crack length which will result in instability at the 
same point on the R-curve, and each combination has an associated fracture stress. To 
cause a small change (dA c ) in the instability point, there must be a change in the speci- 
men width or the initial crack length or both, and there will be a change in the fracture 
stress. Thus, when differentiating equation (2), the terms W, a, and a Q are treated 
as variables. First, note that 


dY _ _9Y_ + _3Y dA , 8Y dW 

dA 3a dA 3A dA 3W dA 
o 


n_ dY 
W dX 



1 - 


a + A .IT 

o c dW 
W dA, 


Now, differentiating both sides of equation (2) with respect to A c yields 


— E'G = Y 2 a 2 2 a /—£■+ 1 - ^2 — _dW 


dA 


c c 


: V dA c 


W 


dA, 


! + Y c (a o + A c> — 
c ° c dA 


-fa 2 ) 

V c ) 


and substituting equation (3) results in 


+ 


yV 

c c 



= g' (A c ) 



(i + 2 eg 


da 
o_ 

^c 


2 a 

c 


a + A 
o c 


W 


dW 

dA c 


+ < a o +A c> 



( 5 ) 


At this point it is helpful to reduce the problem to one of a single independent variableby 
prescribing the manner in which a Q and W may vary. Three cases will be considered. 


Case I: W = Constant (dW = 0) 


Assume that there is a function f such that we can define 
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f(a ) = c 
' o c 




f'(a 


,- d (°c) 
> da. 


> for W = constant 


Then equation (5) becomes 


0 ' [ (1 + 2 “c >f(a 0> + (a o + A c> f, < a o>] 


and, since da /dA / 0, we have 
o c 


0 = (1 + 2a c )f(a o ) + (a Q + A c )f'(a Q ) 


(6a) 


Case II: a Q /W = Constant JdW = (W/a^dal 
Now assume that there is a function h such that 

h(a o ) = <% 


A 


ate) 

h'(a.)= 

da_ 


> for — = constant 
W 


J 


Now equation (5) becomes 


da 

0 = — — 
dA 


1 - 2a„ )h(a o ) + (a D + A c )h'(a o ) 


and since da /dA / 0, we have 
o c 


0 = 


A 

\l - 2a n — )h(a Q ) + (a Q + A c )h'(a o ) 


(6b) 
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Case III: a Q = Constant (da Q = 0) 


This ease is of limited usefulness but is included for completeness. Assume that 
there is a function j such that 


j(W) = a 


A 


j'(W) s 


dW 


>for 


a = constant 
o 


Now equation (5) becomes 


0 = — (a o + A ) 

dA ° 
c 


and since dW/dA c ^ 0 and (a Q + A c ) ^ 0, 


' 2a 

j 

W 


(W) + j'(W) 


2 a 

0 = £ j(W) + j’(W) 

W 


(6c) 


For cracks in infinite bodies, 0 = 0 and equation (6a) becomes 


A 

c 


f < a o> 

-f(a o ) 


a 


o 


(7) 


which, after substituting the infinite-body formulations for f( a Q ) and f' (a Q ) from the 

Appendix, gives A for any a in terms of the empirical parameters. Then terms 
c o 

can be rearranged to give a as a function of A , say 

O V 

a o = F < A c> 


Substituting this function into equation (2) yields 

E' g o =y 2[ F (A o ) + Ac] • 

Since equation (7) gives A for any value of a , equation (8a) must give E’G for any 

CO c 

and all values of A c> which is a definition of the R-curve. Thus, after writing the func- 
tion F in terms of the empirical parameters, it is appropriate to write equation (8a) in 
the general terms of E'G^ and A, rather than E'G c and A c> The end result is an 
explicit ERC formulation in terms of the empirical parameters. 
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To determine the EEC for cracks in finite bodies, an indirect method is required. 
First, the finite-body formulation for f, h, or j and its derivative (see the Appendix) 
are substituted into the appropriate one of equations (6). Because of the more compli- 
cated nature of the finite-body formulations, it is unlikely that an explicit function F(A c ) 
will be obtainable. But for any given value of A c , a Q is a root of equation (6) which 
may be found by standard numerical methods and which represents a single value of 
F(A ). For Case HI, of course, one determines W = F(A ) and solves for W as a root 
of equation (6c). Substituting the root into equation (8a) (if Case I) or into 


E'G = Y 
c c 


[ F «V + < 

• h [ r &v] 

[ F <\P + A J 

• i [ F <\,>] 


if Case II 
if Case III 


(8b) 

(8c) 


as appropriate yields a discrete point on the EEC. By incrementing A £ and repeating 
the calculation, the EEC can be determined point by point. 


ANALYTICAL COMPAEISONS 


Dimensionless Equivalent E-Curves 


It is helpful and more efficient to first compare semiempirical fracture analyses 

(SEFA) and equivalent E-curves (EEC) on an analytical basis. This is most easily done 

using the problem of a crack in an infinite plate as a baseline. 

The infinite-body EEC for Kuhn's analysis (ref. 1) is obtained using the method 

described in the paragraph containing equations (7) and (8a). Substituting equations (Al) 

2 

and (A2) from the Appendix into equations (7) and (8a) with - tt yields 


A 

c 



(9a) 


7TCT^ A 

E'G R y (10a) 

JV o 

1 + C z A 
m 

_l/2 

where C is an empirical parameter having units (L ). Equation (10a) is plotted 
m 

in dimensionless form in figure 5(a). This curve obviously resembles an E-curve and 
might be expected to closely match some (but not all) experimental E-curves. 
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The infinite-body ERC for Orange's analysis (ref. 2) is obtained in the same man- 
ner. Using equations (A3) and (A4) results in 


A =• 
c 


J R 


K \ 2 



01. / 

i u/ 

(9b) 

7T 



(10b) 


where K u is an empirical parameter having units (FL -3 / 2 ). These equations define a 
single point. In order to relate this single point to the R-curve concept, the point may 
be thought of as the corner of a step-function, and that step-function might in turn be 
considered as a very simple approximation of an actual R-curve. 

The infinite-body ERC for Feddersen's analysis (ref. 3) is obtained by substituting 
equations (A5) and (A6) into equations (7) and (8a) as before, yielding 



for a Q < (9/47r)(K c /<7y g )^ where K c is an empirical parameter having units (FL -3//2 ). 
Equation (9c) requires that critical crack extension decrease as original crack length in- 
creases from zero. This is in direct opposition to the R-curve concept and is not sup- 
ported by any data known to this author. Equation (10c), which is plotted in figure 5(b), 
does not look at all like an R-curve but does satisfy the requirements of coincidence and 
tangency. For a Q = 0, the point of tangency is the right terminus of the curve. As a Q 
increases, the point of tangency moves downward and leftward along the curve. Finally, 
at a Q = (9/47r)(K c /cTy S ) 2 , the point of tangency is the left terminus. 

The infinite-body ERC for Newman's analysis (ref. 4) is obtained by substituting 
equations (All) and (A12) into equations (7) and (8a) as before, yielding 


A 


c 



(9d) 
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2 


(10 d) 


E ' G R 


o 2 A 
u 


m 

7T 


K* 


-3/2 

where is an empirical parameter having units (FL ' ) and m is a dimensionless 

empirical coefficient which is not greater than unity. Equation (lOd) is plotted in di- 

o 

mensionless form in figure 5(c). This ERC is asymptotic to E'G-^ = K^, and the coef- 
ficient m determine^ the rapidity of the approach. As m decreases from unity to 
near-zero, the ERC develops a progressively sharper knee. The flexibility of this 
two-parameter ERC should allow it to match R-curves for a wide range of real ma- 
terials . 

The infinite-body ERC for Bockrath's analysis (ref. 5) is obtained by substituting 
equations (A15) and (A16) into equations (7) and (8a) as before, resulting in 


A 

c 



(9e) 


cu/(2+o;) 

E ' G R=f < 2 + “> K Tc(| A ) < 10e > 

where w is a dimensionless empirical coefficient and Krp c is an empirical parameter 
having irrational units (FL^). Equation (lOe) is plotted in figure 5(d). This ERC has 
no asymptote, and its slope is infinite at A = 0. Except for notation, it is identical to 
the R-curve model proposed by Broek as equation (10) of reference 9. Broek's model 
was derived using R-curve concepts and the experimental observation that, for small 
cracks in wide specimens, the critical crack length is often proportional to the initial 
crack length. 

At this point we can state the following. For each SEFA, in its infinite-body form 
at least, there is indeed an ERC. For four of the analyses considered, the ERC re- 
sembles or approximates an actual R-curve. The ERC for Feddersen’s analysis does 
not resemble an R-curve and will not be considered further. 


Comparisons Using Synthetic Data 

Hypothesis (2) of the INTRODUCTION postulates that a SEFA will correlate residual 
strength data closely if its ERC closely matches the actual R-curve and will correlate 
poorly if the match is poor. To test this hypothesis we would need, as a minimum, both 
residual strength data (over a wide range of crack lengths) and R-curve data for two 
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materials having significantly different R-curve shapes . Further, it would be desirable 
to have data for several specimen sizes (including quasi-infinite) and specimen types. 
Since no such body of data is known to this author, it was necessary to synthesize one. 
This was done by formulating two R-curve equations using the following guidelines . 

First, to avoid exact fits, neither equation should be mathematically equivalent to one of 
of the ERC formulations previously derived. Second, one curve should have a definite 
knee, the other should be gently curving. Using these two R-curve equations, synthetic 
test data can be calculated by instability analysis for any specimen size and type. A 
practical advantage of this approach is the total absence of data scatter. 

Unobtainium is assumed to be a heat- treatable material. In the annealed condition, 
its ultimate tensile strength is 150 and its R-curve is given by 

E’G R = 8000 VlOA - A 2 (11a) 

Since the material is imaginary, the units will be left to the reader’s imagination. In 
the aged condition, its ultimate tensile strength is 200 and its R-curve is given by 

E'G d = 50 Q P° arctan (10 A) (lib) 

R 7 r 

These are shown in figure 6. The coefficients in equations (11) were selected so that 
the significant features of both curves would lie within the ranges 0 ^ E'G R < 25 000 
and 0 < A < 1. 

The pseudotest data points are calculated using conventional instability analysis as 

described in the paragraph containing equation (4). The ’’specimens” that are studied 

here were sized as follows. For the infinite-width pseudotests, the initial crack half- 

lengths a were chosen (by trial and error) to give A values well distributed over 
o c 

the entire R-curve. For the finite-width center-crack pseudotests, the same initial 
crack lengths were used but the specimen widths were fixed at eight times (first ser- 
ies) and four times (second series) the largest initial crack half-length. For the com- 
pact specimens, the ratio a Q /W was fixed at 0.5 and the widths chosen (again by trial 
and error) to give approximately the same A c values as those calculated for the 
infinite-width pseudotests. The calculated values of stress and crack extension at in- 
stability are given in table I. The values of stress and initial crack length were then 
used as inputs to the various semiempirical analyses. 

The empirical parameters were determined as follows. Kuhn's parameter C m 
was calculated for each specimen in the infinite-width series using equation (Al) . The 
simple average of seven values C m is given in table I. The bar is used here to denote 
the average value for one data set. Orange's parameter was also calculated for 
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each specimen in the infinite-width series using equation (A3). The average value 
given in table I is a weighted average determined in the same manner as equation (6) of 
reference 2. Newman's parameters and in were determined using the least- 
squares procedure given in Appendix C of reference 4. Bockrath's parameters K 
and co were, determined by a least-squares fit of equation (A15). Since Bockrath's 
method is restricted to cases where the crack area is less than 10 percent of the gross 
area, specimens having a Q > W/20 were excluded from the fit. 

The equivalent R-curves were calculated as follows. For the infinite- width series, 
the empirical parameters from table I were simply substituted into the appropriate one 
of equations (10). For the remaining series, the indirect method (described in the par- 
agraph containing eqs. (8b) and (8c)) was used. Specifically, for the finite-width center- 
crack series, equations (A9) and (A10) or equations (A15) and (A16) were used along 
with equations (6a) and (8a) and the appropriate empirical parameters. For the compact 
specimen series, equations (A13), (A14), (6b), and (8b) were used. 

When determining the Newman ERC for constant- finite-width specimen configura- 
tions, some precautions must be taken. For values of A c near zero, there are three 
positive roots of equation (6a); one at A Q near zero, one at A q ~ 0.44, and one at A 
near unity. The middle root gives an ERC which resembles the Feddersen ERC of fig- 
ure 5(b) and which should be ignored. The other two roots give ERCs which are numer- 
ically similar but distinctly different. As discussed earlier in the section R-CURVE 
CONCEPTS, the R-curve concept implies a particular relation between fracture stresses 
for short cracks and long cracks. Since Newman's analysis does not match this relat- 
tion, we have separate ERCs for short and long cracks. As the value of A c is in- 
creased, the middle root converges upon one of the other roots whereupon the converged 
roots vanish. In this report, the Newman ERC for constant-finite-width configurations 
is determined using the least positive root of equation (6a) over the range of A c /W for 
which three positive roots exist. 

The residual strength of the infinite-width Unobtainium is shown in figure 7. For 
the annealed condition, Bockrath's semiempirical fracture analysis (SEFA) provides a 
nearly perfect fit to the pseudodata. When ranked according to the sum of the squares 
of the deviations, Newman's SEFA, Kuhn's, and Orange's follow in that order. For the 
aged condition, the ranking is quite different. Here Newman's SEFA provides a nearly 
perfect fit, with Orange's, Kuhn's, and Bockrath's following in that order. The equiva- 
lent R-curves (ERCs) are shown in figure 8. For the annealed condition, the Bockrath 
ERC matches the actual R-curve almost perfectly. When ranked according to the inte- 
gral of the square of the deviation, the Newman ERC, the Kuhn ERC, and the Orange 
ERC follow in that order. For the aged condition, the Newman ERC is the best match 
to the actual R-curve. The Kuhn ERC, the Bockrath ERC, and the Orange ERC follow 
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in that order. The Orange ERC, although crude, is a better match to the aged mate- 
rial than to the annealed. 

Critical crack extensions derived from the Newman and Bockrath ERCs are shown 
in figure 9. The Bockrath ERC predicts crack extension fairly closely for the annealed 
condition, but neither ERC gives a very good prediction for the aged condition. This 
last observation should not be surprising. Returning to figure 2, it should be apparent 
that a small change in the slope of the R-curve near the instability point will cause a 
large change in the critical crack extension but only a small change in the critical 
stress. Thus when instability calculations are made from any R-curve, be it actual or 
equivalent, fracture stress can be predicted with much more confidence than can crack 
extension. 

Since the Orange ERC is rather crude and since Kuhn's SEFA is equivalent to a 
special case of Newman's (see the Appendix, following eq. (A12)), these two were not 
considered further. The Newman and Bockrath ERCs for the finite-width center-crack 
series are shown in figures 10 and 11. Here the same trends are seen as in the infinite- 
width series. The Bockrath ERC is the better match for the annealed condition, while 
the Newman ERC is the better match for the aged condition. Note in table I that the em- 
pirical parameters Kr rc , a?, Kj., and m all vary slightly with specimen width. The 
ERCs shown in figures 8, 10, and 11 are also distinctly different for different specimen 
widths, but the differences are slight. 

Since Bockrath's SEFA does not include the compact specimen geometry, only the 
Newman ERC was determined for the compact specimen series. Note again in table I 
that Newman's empirical parameters for the compact specimen series differ from those 
obtained for the various center-crack series. The differences in the parameter 
are slight; the differences in m are somewhat larger. The Newman ERCs for the 
compact specimen series are shown in figure 12. The correspondence between equiva- 
lent and actual R-curves is not quite as close as for the center-crack series, but the 
general trend is the same. A much better match is obtained for the aged condition than 
for the annealed. 

The results of this exercise using synthetic data can be summarized as follows. 
Hypothesis (2) of the INTRODUCTION postulates that an SEFA will correlate residual 
strength data closely if its ERC closely matches the actual R-curve and will correlate 
poorly if the match is poor. Strictly speaking, this hypothesis cannot be proven, since 
the ERC magnitude and shape depend on empirical parameters which must be obtained 
from residual strength data. However, the converse appears to be true. That is, if an 
SEFA correlates residual strength data closely, its ERC will closely match the actual 
R-curve. Furthermore, it is apparent that if, for a given material and thickness, the 
R-curve is unique the various empirical parameters are not, and vice versa. Also, al- 
though it is possible to predict the critical crack length in terms of the original crack 
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length and the empirical parameters, this cannot be recommended. The crack extension 
at instability is quite sensitive to small differences in the slope and magnitude of the 
R-curve near the instability point. 


COMPARISONS USING ACTUAL TEST DATA 

As mentioned earlier, experimental studies containing both residual strength and 
actual R-curve data are relatively few in number. Nevertheless, enough were found in 
the literature to allow some comparisons to be made using actual data obtained from 
real materials. 


NASA Data for 2014-T6 Alu m inum Alloy 


In reference 10 this author presented test data for 2014-T6 aluminum alloy speci- 
mens 1.5 millimeters (0.06 in.) thick, tested at 77 K (-320° F). Figure 14 in refer- 
ence 10 presented typical curves of crack growth against applied stress for notches of 
six initial lengths in 30-centimeter (12-in.) wide specimens. Those curves were de- 
veloped by plotting individual crack growth data points for replicate specimens, then 
drawing a smooth curve to give a good visual average. For the present report, those 
data were reanalyzed. The crack extension resistance and effective crack length were 
computed for each data point as 


E'Gr = o 2 


r cififf secant 



a eff “ 



respectively. Since these equations are transcendental, an iterative solution was re- 
quired. A total of 176 data points were obtained from 17 specimens with initial crack 
lengths 2a ranging from 3 to 100 millimeters (l/8 to 4 in.). The empirical parame- 
ters for the Newman and Bockrath SEFAs were determined in the manner described 
earlier. As before, only specimens with a Q < W/20 were included in the Bockrath 
analysis. The fitted empirical parameters are listed in table II. 

Residual strength is shown in figure 13. Newman's SEFA gives a good fit over the 
entire range. Bockrath's SEFA fits the short-crack data fairly well, but the fit would 
be poor if extrapolated to longer cracks. The R-curve data points and the equivalent 
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R-curves are shown in figure 14. Both ERCs fit the data rather well, with Bockrath's 
somewhat better at small crack extensions and Newman's somewhat better for larger 
extensions. These results suggest that the ERC concept applies to real data as well as 
to synthe sized data. 


Battelle Data for 707 5- T7 351 Aluminum Alloy 

Reference 11 contains test data for 7075-T7351 aluminum alloy specimens 25 milli- 
meters (1.0 in.) thick and 40 centimeters (16 in.) wide. This data set is unusual in that 
the ratios of initial crack length to plate width range uniformly from near-zero to near- 
unity; also, the specimens were instrumented to measure crack opening displacement 
(COD) during loading, and the stress corresponding to a 5-percent secant offset on the 
load-COD curve is reported. A secant offset corresponds to an increase in the effective 
crack length. From an R-curve one should be able to predict the stress required to 
cause such an increase. 

Since Bockrath's analysis is restricted to small cracks, only the Newman ERC was 
used here. First Newman's empirical parameters were determined by fitting the re- 
sidual strength data as described earlier. The stress at 5 percent secant offset is cal- 
culated as follows. According to reference 12, the crack opening displacement 6 for a 
finite-width center-crack plate is 


6 = 





where 


V^X) = -0.071 - 0.535 X + 0.169 X 2 + 0.020 X 3 - 1.071 A -1 ln(l - X) 
The slope of the load-COD curve is thus proportional to 



Let A 5 be the crack extension at 5 percent secant offset. Then a crack of half- 
length a Q + A 5 will give a slope 5 percent less than a crack of half-length a Q , or 

h + a 5 > v i < v ]’ 1 = °- 95 [ Viv ]' 1 
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where X = 2 (a + AA/W. Rearranging terms results in 
o o o 

, V (X ) a + A k 
0 = 1 1 o' _ o 5 

0.95 V 1 (\ 5 ) a Q 

and A r is the least positive root of the previous equation, which may be found by stan- 
dard numerical methods . Now that A,- is known, the corresponding crack extension 
resistance E’Gg is obtained from the Newman ERC. Finally the stress at 5 percent 
secant offset cr_ is 


X E'G /?r(a + A_) secant 7r(a + AA/W 

a 5 = I § 2 2 ° 13 (12) 

1 - 0.025 Ag + 0.06 A.g 

Here the modified-secant stress intensity factor expression from reference 12 was used 
for greater accuracy at relative crack lengths near unity. 

The fitted empirical parameters are given in table II. Residual strength is shown 
in figure 15, where Newman’s SEFA can be seen to fit the data reasonably well. The 
stress at 5 percent secant offset is shown in figure 16. Except for the shortest crack, 
the stress predicted by equation (12) agrees well with the reported data, the scatter be- 
ing no worse than the scatter in the residual strength. Since the length of the shortest 
crack is less than twice the plate thickness, it is likely that the two-dimensional repre- 
sentation implicit in the analyses is no longer valid. The otherwise good agreement 
suggests that, even though the actual R-curve is not available for comparison, the New- 
man ERC must closely approximate it. 


Alcoa Data for Compact Specimens 

Reference 13 contains test data for compact specimens of 2219-T851 aluminum alloy 
in four thicknesses. Specimen sizes were varied, but the initial crack length was al- 
ways half the specimen width. The maximum loads and the loads designated Pq are 
reported. The load Pq is the load corresponding to a 5-percent secant offset on the 
load-COD curve and is used to compute the plane strain fracture toughness K^. As 
was done for the Battelle center-crack specimens, one should be able to predict these 
loads from an ERC. 

Since Bockrath's SEFA has not been formulated for the compact specimen geometry, 
only Newman’s SEFA is examined here. Newman's empirical parameters were deter- 
mined for each thickness in the manner described earlier for compact specimens. 
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According to reference 14, the crack opening displacement V Q measured in a Kj c 
test is given by 


2E'V 

2 = 120.7 

P 


1065.3 X + 4098.2 7? - 6688.0 X 3 + 4450.5 A 4 


Using this expression, the crack extension corresponding to a 5-percent secant offset 

with a/W = 0.5 is A,. = 0.023 a and the corresponding crack extension resistance 
c o o 

E'Gg is thus obtainable from the Newman EEC for the appropriate thickness. The load 
Pq is then 



(13) 


where B is the specimen thickness and Y^ is the stress intensity calibration factor 
corresponding to a = a Q + = 1.023 a Q . The equation from reference 8 given follow- 

ing equation (4) in this report gives Y^ = 13.79 for X = 0.5115. 

The predicted loads Pq for each thickness are shown in figure 17 . Agreement 
with the reported test loads is generally good, with the average error for the 32 tests 
being about 8 percent and the largest absolute error being less than 9 KN (2 kip). This 
figure indicates that the Newman ERC can also approximate the actual R-curve using 
residual strength data from compact specimens. 


Boeing Data for 2219-T87 Aluminum Alloy 

Earlier it was shown using synthetic data that, if the actual R-curve is unique, one 
obtains slightly different values of the empirical parameters from data sets for speci- 
mens having different widths. Data from the Boeing Co. for 2219-T87 aluminum alloy 
specimens have the same characteristics. These data originally appeared in an internal 
report (Eichenberger, T. W.: Fracture Resistance Data Summary. Report D2-20947, 
Boeing Airplane Co., June 1962) but are also tabulated in reference 1. Center-crack 
specimens 2.5 millimeters (0.10 in.) thick were tested. The data for specimens 60 and 
120 centimeters (24 and 48 in. ) wide are used here because they cover a wide range of 
crack lengths. 

Bockrath's analysis was not applied since only one of the wider specimens had 
a Q < W/20 . Newman's parameters were determined separately for each specimen width, 
and somewhat different values were obtained as can be seen in table II. The residual 
strength curves fit the data quite well, as can be seen in figure 18, with the average 
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error being less than 3g percent. Using a method that is outside the scope of this re- 
port, it was found that the actual R-curve for this material could be estimated by 


E'G r = 8.07X10 15 A 0 ' 554 

o q 

where E'G R is in N /m and A is in cm, or by 

E'G r = 11.2X10 9 A 0 * 554 

2 S 

where E'G„ is in lb /in. and A is in inches. Residual strengths calculated from 

XV 

this equation using conventional instability analysis are also shown in figure 18 . The 
agreement is very slightly better than for Newman's SEFA, the average error being less 
than 3 percent. The estimated R-curve and the Newman ERCs are shown in figure 19, 
and as expected the differences are small. Although these data tend to support the con- 
cept of a unique R-curve, the differences are so small as to be within the bounds of 
probable data scatter. 


CONCLUDING REMARKS 

The relations between several semiempirical fracture analyses (SEFA) and the 
R-curve concept of fracture mechanics are examined and the conditions for equivalence 
between a SEFA and an R-curve are derived. A hypothetical material is employed to 
study the relation analytically. Equivalent R-curves are developed for several real ma- 
terials using data from the literature. The results of this study lead to the following 
conclusions: 

1. For each semiempirical fracture analysis (SEFA) there is an equivalent R-curve 
(ERC) whose magnitude and shape are determined by the SEFA formulation and its em- 
pirical parameters. The ERC is equivalent in that it predicts exactly the same relation 
between fracture stress and initial crack length (residual strength) as the SEFA. 

2. If, for a given set of data, a SEFA correlates residual strength closely, its ERC 
will closely approximate the effective R-curve of the material. 

3. It is possible to predict the critical crack length in terms of the initial crack 
length and the fitted empirical parameters. However, this is not recommended since 
the crack extension at instability is particularly sensitive to the slope of the R-curve. 

4. Of the five SEFAs examined, Newman's appears to be the most generally useful. 
Bockrath's SEFA, which is only formulated for quasi-infinite bodies, is too restrictive 
for widespread use. Three other SEFAs do not appear to warrant further consideration. 

5. If the effective R-curve is unique, then the various empirical parameters cannot 
be constant, and vice versa. 

20 
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The analytical comparisons made herein indicate that the variations in Newman’s 
parameters are small enough that the differences may well be within the range of nor- 
mal data scatter for real materials. Thus a very carefully planned and conducted ex- 
periment would be required to determine which concept (R-curve or SEFA) is more uni- 
versally applicable. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 18, 1979, 

505-02. 
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APPENDIX - SEMIEMPIRIC AL FRACTURE ANALYSES 


In this Appendix, equations from references 1 to 5 are rewritten using the notation 
of this report. Note the following definitions from the text: 


« a o> = 4 


f, <v $ 




> for W = constant 


da. 






Ma o ) - og 


h’(aj 


y for — = constant 
W 


da. 


P. Kuhn (1968) 

Equations (3) and (4) of reference 1 give the fracture stress for a finite- width 
center-crack plate as 


a c = u 


u * 1 - V [i + C m - V * 1 + VJ 


-1 


_l/2 

where C m is an empirical parameter having units (L ' ). Thus for an infinite plate, 


f(a o ) 


2 / 
7 1 

U \ 


1 + c. 


m 



f, ( a 0 ) = ' f(a o ) 


1+ C 


m 



(Al) 


(A2) 
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T. Orange (1969) 


Equation (8) of reference 2 gives the fracture stress for a finite-width center-crack 
plate as 


a 



W tan 



where is an empirical fracture toughness parameter having units (FL ) . For an 
infinite plate, this reduces to 


f(a ) = K 
' o' u 


7ra + 

o 



(AS) 


f '(a o )=- f ( a o ) 



-1 


(A4) 


C. Feddersen (1970) 

In reference 3, the fracture stress for a finite-width center-crack plate is given in 
equations (6) and (10) by 


a = o 
c ys 



for a 

o 



and in equation (7) by 


a = K (7ra ) 1//2 
c c ' o' 


for 



W 

6 


and in equations (8) and (11) by 
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°c=f tL -V 



for a 

o 


W 

6 


where a is the material's yield strength and K is an empirical fracture toughness 

yS C 

parameter. For an infinite plate, the first of these three equations reduces to 


f(a o ) 


°ys 



f'(a o )-- f ( a o) 




and the second to 




f’(a o ) = 


- f (a o ) 


a 

o 


(A5) 


(A 6) 


(A7) 

(AS) 


J. Newman (1972) 

Equation (12) of reference 4 gives the fracture stress for a finite-width center- 
crack plate as 



where is an empirical fracture toughness parameter and m is a dimensionless 
empirical coefficient which is not greater than unity. Rewritten and differentiated, this 
is 
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_ -2 



Note that, if we let m = 1 and k f = cr u \n/ c m > equation (All) reduces to equation (Al). 
For the compact (tension) specimen, equations (1), (7), and (16) to (19) of reference 15 
result in 


i-2 


h(a Q ) = K| 


mK, 


a + 
o 


L(i-V u u 


(A13) 


h '(a o ) 


- h < a G ) 


o y o 


l a 


V 



(A14) 


where 
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L = 



G. Bockrath (1972) 

Equation (1.3) of reference 5 gives the fracture stress for a center-crack plate as 


'Tc'V 

where w is a dimensionless empirical coefficient and K-p c is an empirical toughness 
parameter having irrational units of (FL 60 ). This equation is limited to A. £ 0. 1, which 
approximates an infinite plate. Thus, 

f( a o ) = 4 c (a o r 2/(2+a,) (A15) 


f'K) = - 


f(a ) 
o 


2 + co a 


(A.16) 
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TABLE I. - PSEUDOTEST DATA FOR HYPOTHETICAL MATERIAL UNOBTAINIUM 


AND FITTED EMPIRICAL PARAMETERS 


(a) Annealed condition; ultimate tensile strength cr^ = 150 


K 


K, 


Tc 


K, 


a o 


Center crack specimens 


Compact specimens 


W = 

1 

= OO 

W = 

8.8 

W = 

4.4 

( a 0 = 0 . 5 W) 


A 

(7 

A 

cr 

A 

(7 

W 

A 



c 

C 

c 

c 

c 

C 


c 

c 

0.10 

0.0980 

112.56 

0.0971 

112.42 

0.0945 

112.00 

1.5 

0.1031 

6.073 

.21 

.2015 

93.25 

.1939 

92.75 

.1757 

91. 36 

3.0 

.2044 

5.094 

.32 

. 3008 

83.70 

.2773 

82.71 

.2313 

80.08 

4.5 

.3037 

4.591 

.44 

.4044 

77.08 

. 3529 

75.43 

.2692 

71.29 

6.0 

.4012 

4.261 

.56 

.5036 

72. 36 

.4137 

69.98 

.2901 

64.26 

7.5 

.4969 

4.019 

.80 

.6897 

65.82 

.4996 

61.79 

.2986 

52.84 

11.0 

.7126 

3.630 

1.10 

.9016 

60.36 

.5580 

54.09 

.2713 

41.02 

14.0 

.8893 

3.400 

C 

1. 

355 








m 
















inol computea 




K u 

111.8 








Etc 

62.13 

60 

.76 

59 

.52 




U) 

1. 

853 

1 . 

727 

1 . 

642 




K f 

165.7 

172.8 

159.6 


157.7 


m 

0.8549 

0 . 8497 

0.7800 


1.0 


(b) Aged condition; 

ultimate tensile strength 

c u = 200 


a 

o 


Center crack specimens 


Compact specimens 


0.05 

.36 

1.0 

2.0 

3.2 

6.4 

12.0 


W = 


c 

0.0953 
.2014 
. 3067 
.4135 
.5093 
.6978 
.9346 


162.93 
100.08 
69.77 
52.91 
43. 37 
31.93 
23.95 


1.820 


140.7 

63.80 

0.838 

163.3 

0.7266 


W = 96 


c 

0.0953 

.2014 

.3064 

.4121 

.5054 

.6791 

.8566 


c 

162.93 

100.07 

69.74 

52.82 

43.21 

31.50 

22.87 


W = 48 


c 

0.0953 
. 2013 
. 3055 
.4081 
.4942 
.6286 
.6760 


c 

162.92 
100.05 
69. 64 
52. 58 
42.73 
30.21 
19.60 


W 


1.2 

4.0 

8.0 

15.0 

23.0 

45.0 

75.0 


0.1011 

.2049 

.2943 

.4042 

.5001 

.6971 

.8971 


8.006 

5.864 

4.592 

3.579 

2.988 

2.221 

1.759 


Not computed 


66.57 

1.151 

166.0 

0.7378 


68.38 

1.310 

168.4 

0.7469 


167.0 

0.9533 
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TABLE II. - EMPIRICAL PARAMETERS FOR TEST DATA. FROM THE LITERATURE 


Data 

source 

Alloy 

Specimen 

constant 

Empirical parameters 

(ref. 4) 

Empirical parameters (ref 

K f 

m 

K Tc 




MN-m"^ 2 

ksi V in. 


MN-m 1 ' 761 

lb-in.- 1 ' 761 

Ref. 10 

2014-T6 

W = 30 cm 

89.9 

81.8 

0.8094 

99.70 

34 790 

Ref. 11 

7075- T7 351 

W = 40 cm 

84.1 

76.5 

1.0 



Ref. 13 

22 19- T8 51 

B = 13 mm 

101.4 

92.3 

1.0 





B = 25 mm 

78.2 

71.1 

1.0 





B = 38 mm 

65.4 

59.5 

.8921 

Not computed 



B = 50 mm 

58.4 

53.1 

. 6997 



Ref. 1 

2219-T87 

W = 60 cm 

183.1 

166.6 

0.8841 





W = 120 cm 

210.2 

191.3 

1.0 




a Only specimens having a Q £ W/20 were included in least squares fit. 



CRACK LENGTH, a 


Figure 1. - Schematic representation of R-curve instability concept 
for an infinite body. 
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Figure 3. - R-curve instability for a wide range of initial crack lengths. 
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Figure 4. - Dimensionless residual strength for center-crack specimens 
of a material with an R-curve of the general form E'G R ■ A(A) b . 
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Figure 5. - Dimensionless R-curves equivalent to various semiempirical fracture analyses for the case of a crack in an infinite plate. 
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(a) Annealed condition: o u =150: 


(b) Aged condition: o u =2D0. 


Figure 7. - Residual strength of hypothetical material Unobtainium, infinite width series; various semiempirical fracture 
analysis (SEFA) fit to pseudotest data. 
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CRACK EXTENSION RESISTANCE, E'G 




(b) Aged condition. 

Figure 8. - Actual R-curves and equivalent R-curves 
(ERC) for hypothetical material Unobtainium; infinite 
width series. 


CRACK EXTENSION, 



(a) Annealed condition. (b) Aged condition. 

Figure 9. - Critical crack extension for hypothetical material Unobtainium; infinite width series; calculated from actual R-curves and from equiva 
lent R-curves (ERC). 
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(a) Width W = 8. 8. 
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CRACK EXTENSION, A 


NEWMAN ERC 
BOCKRATH ERC 
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(b) Width W = 4. 4. 


Figure 10. - Actual R-curve and equivalent R-curves (ERC) for hypothetical material Unobtainium ; finite width; annealed 
condition. 
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Figure 15. - Residual strength of 7075-T7351 aluminum alloy plate (ref. 11)-, 
Newman's semiempirical fracture analysis (SEFA) fit to test data. 
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Figure 16. - Stress at 5 percent secant offset 05 for data from reference 11, 
predicted using equation (12). 
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Figure 18. - Residual strength of 2219-T87 aluminum alloy sheet 
(ref. II calculated from Newman's semiempirical fracture analy- 
sis (SEFA) and from estimated R-curve. 
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Figure 19. - Equivalent R-curves (ERCI and estimated R-curve for 
2219-T87 aluminum alloy sheet (ref. 1). 
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